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Abstract 

The dynamic structure factor S'(q, u) of a harmonically trapped 
Bose gas has been calculated well above the Bose-Einstein con- 
densation temperature by treating the gas cloud as a canonical 
ensemble of noninteracting classical particles. The static struc- 
ture factor is found to vanish oc in the long- wavelength limit. 
We also incorporate a relaxation mechanism phenomenologically 
by including a stochastic friction force to study ^(q, a;). A signif- 
icant temperature dependence of the density-fluctuation spectra 
is found. The Debye-Waller factor has been calculated for the 
trapped thermal cloud as function of q and the number J\f of 
atoms. A substantial difference is found for small- and large-A/" 
clouds. 
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1 Introduction 



There has been much interest in the study of trapped atomic gases during 
more than a decade. [21 E] Experimentahsts have observed not only density 
profiles of these systems but they have successfully detected particle-number 
and -density fiuctuations in ultra-cold atomic gases. [1] Progress has been 
made in understanding spatial and temporal correlations of density fiuctua- 
tions in a variety of inhomogeneous cold-atom systems O El O [8] for T=0 
as well as below the Bose-Einstein condensation (BEC) temperature Tg. In 
addition there have been extensive studies of dynamical density fiuctuations 
for homogeneous systems. [9l [TO] 

The purpose of this work is to calculate the diagonal elements of the ma- 
trix of time-dependent correlation functions of density fiuctuations for a 
harmonically trapped classical ideal gas, which determine the coherent scat- 
tering properties of the system. Motivation for calculating the dynamical 
structure factor is its expected strong temperature dependence which could 
be of use in estimating the cloud temperature experimentally. Moreover, 
the exactly solvable model not only provides the limiting case of a system 
in the absence of interactions but also furnishes a useful result for trapped 
dilute gases well above the BEC temperature. Justification for a classical 
treatment is based on the fact that for an atom trap of typical trap fre- 
quency f2=120 X 27rs~^ containing A/'=2 x 10^ rubidium atoms which implies 
a critical temperature Tc=^f2A/'^/^/[C(3)^/^A;B] ~ 147 nK of Bose-Einstein 
condensation [2l Chap. 2.2], the ratio HQ over k^T is very small in the 
temperature range well above Tc, where many experiments are performed, 
because ?iVL/{k-QTc) ~ 0.04. Moreover, at these temperatures the reduced 
density p=nX^ ^ 0.01 is small implying A ~ 0.008 L, i.e. the deBroglie wave- 
length A is still very small compared to the linear extension L of the atom 
cloud. 

In addition to the restoring force acting on each particle due to the trap 
potential, we allow for a frictional force when studying the dynamical be- 
havior of the trapped gas. Simultaneously, a time-dependent random force 
is assumed to act on each oscillator. This stochastic force will keep particles 
moving, thus preventing them from being slowed down by friction and com- 
ing to a rest in the trap-potential minimum. Thermal averages are therefore 
understood to include additional averaging over stochastic variables. In this 
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way, the assumption of time-translational invariance will remain justified 
throughout in spite of frictional forces. 

In Sec. m we present the basic definitions and the correlation functions to 
be studied. In Sec. El the one-and two-particle densites and the static 
structure factor are obtained. In Sec. HI the intermediate scattering function 
has been calculated in the presence of damping including a stochastic force 
with Gaussian noise. Resulting low-friction density-relaxation spectra are 
discussed as function of temperature. In Sec. [5l density excitations for 
the limiting case of zero friction and the Debye-Waller factor (DWF) are 
discussed. Finally, concluding remarks are given in Sec. O 



2 Theoretical Considerations 

In the description of dynamical properties of a many-particle system, the 
correlation function of density fluctuations 6N{r,t), 

F(r,r',t) = -l(5iV(r,t)5iV(r',0)) , (1) 

and its spatial Fourier transform 

Fqq,(t) = Jd\ Jd\' e-'^-^+'^'-^'F{r,r',t) = ^{6N^{t)6N^i0r) (2) 

play an important role. Here (...) denotes a thermal equilibrium average, 
and the dynamical variable A^(r, t)= Z^jli ^(r ~ ^j{t)) describes the particle 
density at space point r and time t, which, when integrated over all space, 
sums up to the total number of particles, N'= J d^r A^(r, t). Correspondingly, 
the density variable in Fourier space is given by 

AT 

N^{t) = / dV e-^'i" N{r, t) = J2 e-'^i'^^*) (3) 

i=i 

with fluctuation 6Nci{t). The thermal fluctuation of a dynamical variable 
A(t) is defined as 6A{t) = A{t) - {A{t))=A{t) - {A), where invariance 
of the thermal equilibrium state with respect to time translations renders 
{A{t) ) independent of time t. 
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The diagonal elements F(q, t)=Fqq(t) are known as intermediate scattering 
function, because they determine the double-differential coherent-scattering 
cross section observed when scattering neutrons or photons from the system. 
The cross section is proportional to /^dt exp{ituj) {Nq(t)Nq{0)* ). It is 
conventionally written as [TT] 



i9 \ coh 

dV ^ pf 2 



\dujdQj P 
by introducing the van Hove function 



( N W"^ 



(4) 



1 /-oo 

^(q, u:) = — dt e'*'^F(q, t) . (5) 

ZTT J-oo 

In the coherent-scattering cross section Eq.(j4j), /iq=pj — p/ and huj=Ei — Ef 
denote the momentum and the energy, respectively, which is transferred in 
a scattering process from the projectile (neutron, e.g.) to the target system 
(gas), and a denotes the coherent scattering length of the gas. 

We note that the right-hand side of Eq.(jl]) will represent a decomposition 
into elastic and inelastic scattering contributions if, and only if, -F(q, oo)=0 
implying an ergodic many-particle system. For non-ergodic systems with 
F(q, oo) > 0, the van Hove function ^(q, a;) itself will have an elastic peak 
of strength -F(q, oo) in addition to inelastic (a; ^ 0) contributions. As will 
be shown below, the harmonically trapped ideal gas represents not only an 
inhomogeneous but also a non-ergodic many-particle system. 

In general, the elastic coherent-scattering intensity is measured by the DWF 
defined as the ratio of elastic and total coherent-scattering intensity, 



!"« iW) |{iv,>|VAf + F(q,o) ■ 

Here 

/oo 1 / \ 

^d^5(q,^)=5(q) = -^(|5iVqp) (7) 

is the well known static structure factor S{q). We emphasize that no as- 
sumptions regarding translational invariance of the target system have been 
made in the above equations. We note that in Eq.([n]) there are two contri- 
butions to the DWF, in general. It is worth recalling that for a cristalline 
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solid the first contribution is nonzero, while the second, i.e. F(q, cx)), van- 
ishes since the cristalline solid is an ergodic system. On the other hand, for 
a non-ergodic (structural) glass -F(q, oo) > 0, and ( ) is zero for all q 7^ 
in the homogeneous system. In a homogeneous fluid, both contributions are 
zero resulting in a vanishing DWF. However, in the inhomogeneous trapped 
gas studied here, there will be a finite DWF as shown in Sec. below. 

Although Fqq/()f:) for an inhomogeneous system will generally be non-zero if 
q^^q', only the diagonal elements are required for calculating the coherent 
double-differential scattering cross section. That is, the intermediate scat- 
tering function F{q,t)=Ff^q{t) will determine the coherent-scattering cross 
section (aside from the average density, which is responsible for the elastic 
'Bragg-scattering' contribution of strength | {N^) P/A/"). 

We now consider a system of A/" non-interacting particles of mass m moving in 
three-dimensional space under the influence of an isotropic harmonic external 
potential. This collection of oscillators is described by the Hamiltonian 



(8) 



where V{r)=^mri'^r'^ with oscillator frequency Q, and pj and denote mo- 
mentum and position variables of the j-th particle, respectively. For the 
trapped ideal gas described by Eq.([8]), thermal averages of additive single- 
particle variables of the general form A{t) := J2f=i ^{Pj{t),'rj{t)) reduce to 
the simple expression 



( Ait) )=Mld'p T^i-^ / d¥ '"^ ^Jf) A(p(t), f (t)) , (9) 



where pT=mvTi VT=\Jk-BT /m^ and u=vt/^ denote, respectively, the parti- 
cle's thermal momentum, thermal speed, and distance traveled within time 
fl~^. In general, r(t)=R(p, r; t) and p(t)=P(p, r; t) will depend on initial 
momentum p, initial position f , and time t. The vector-valued functions 
R and P will be determined from Hamilton's equations of motion in SecJH 
below. For t=0, however, f(0)=r and p(0)=p. The evaluation of the 6-fold 
integral on the right-hand side of Eq.Q is straight forward then. The static 
averages entering Eqs. (jH O [7]) are discussed in Sec|3l 
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3 Static correlations 



Using Eq.Q, one obtains for the average particle-number density, 

A/'exp I ^ I 

{N{r))= l_'f^ « (iV,)=Are-i''^"\ (10) 



We also note that the mean-square displacement of a particle from the center 
of the trap is 

{'l) = {^pl)=3n\ (11) 

rendering an additional physical interpretation of the length u. According to 
Eqs.f llUHTT]) . in thermal equilibrium the trapped gas will form a spherically 
symmetric cloud in which 90.0(99.0; 99.9)% of the gas particles are contained 
in a sphere of radius a ^2.5(3.4; 4.0)-u around the center of force. 

For the initial value Fqq/(0) describing the static correlations of density fluc- 
tuations, one finds (see also Eqs. (ITBHTTI) below): 

^qq'(O) = ^{SN^SN*^, ) = e-5('?^+^")-' [e"^-^'-' - l) . (12) 

The corresponding result in r-space is given in App|Bl Note the symmetry 
relation, Fqq/(0)=Fq'q(0), and the vanishing small-g limit, limg^o -fqq'(0)=0- 

At this stage, let us point out an interesting fact regarding the static structure 
factor defined in Eq.([7]) which, for the trapped ideal g be read from 

the diagonal elements of Eq. (IT^ : 



S{q) = 1 - e-«'"' . (13) 

As displayed in Fig. [D^a), S{q) q'^u^ + O^q"^) for long wavelengths. Lack- 
ing a well-defined volume occupied by the A/" particles, an overall particle- 
number density n cannot rigourously be defined in the inhomogeneous gas. 
However, the "peak density" 

no:=(A^(r=0)) = ^, L := uV2^ (14) 

presents a reasonable measure of the overall density. Due to the rapidly 
decreasing average density {N{r)) for |r| > u (see Eq.f llOl) ). V=L^ will 



6 



be considered to represent the gas volume. It corresponds to a sphere of 
radius a=3^/'^(7r/2)^/^ l.56u containing about 51% of the total number 
J\f of particles. An estimate of u for ^^Rb in a trap of Q=120 x 27rs~^ at 
T=5Tc ~ 735 nK gives u ^ 11 /im. Assuming the trap to contain A/'=2 x 10^ 
atoms, the peak density will be no ~ 10^^ cm~^ which is typical of the values 
encountered in cold-atom traps. 

For the non-uniform ideal gas considered here, one finds for the average 
two-particle density [TOl Chap. 2.5] 

implying the two-particle distribution function g\/ {t, r') = 1 — 1/Af , which 
is a constant depending on particle number Af. 



4 Dynamic correlations 

It is straight-forward to show that the density-fluctuation correlation func- 
tions of the non-interacting gas take the simple form 

F^.'it) = cil,it)-e'^^^'^'^''^-' (16) 

with Ci'^,(t)=^Ej£i (e-''i-'"^(*)e''i' '"^(°) ) denoting the tagged-particle den- 
sity correlation function, which reduces to 



where {. . .}av explicitly indicates averaging over the stochastic-force variable. 
Inserting into Eq.( |T71) the solution r(t) according to Eq. (]A.3l) with initial time 
to=0 , we obtain 

Csiq,t) = Cj)(t) = exp{-gV[l-0(t)]} (18) 

with initial value Cs{q, 0)=1 implied by Eq. flA.7p . For non-zero friction, the 
long-time limiting value Cs{q, oo) = exp {—q^v?) follows from Eq. flA.Sp . 
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Figure 1: (a) Static structure factor; (b) characteristic frequency for trapped 
gas (full line) and free gas (dashed), and u!=Q (dash-dotted). 



As a consequence, the intermediate scattering function becomes 



F{q,t) = F^^it) 



(19) 



which is displayed in Fig. [2] for various F and one representative wavenumber 
q, will have a damped oscillatory form, in general, with amplitudes ranging 
between 

-Siq)e-^'-' <F{q,t)<S{q) . (20) 

It is readily seen from Eqs.(fT8l).( 1A.8l) that limt^oo F{q,t)=0 for any finite 
friction constant, and oscillations will be over-damped, if F exceeds the 
aperiodic-limit value F=2f2. Evidently, the intermediate scattering func- 
tion of the trapped ideal gas is closely reflecting the motion of a single gas 
particle under the influence of friction and harmonic trap potential. 

We note that the static structure factor S{q) is not affected by the damping 
mechanism introduced above. Further, the short-time behavior. 



F{q,t) = S{q) 



qvT 



(21) 



is governed by Ug, the characteristic frequency defined by the second fre- 
quency moment ujg := Jduj uj'^(f)"{q,uj) of the normalized dynamic structure 
factor, 

0"(^,-) = -%^. (22) 



S{q) 



The static structure factor, which is vanishing as for g — 0, is respon- 
sible for the gap of size Q in the graph of Uq, see Fig. [D^b). Physically, 



8 



F(q,t) 
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Figure 2: Intermediate scattering functions at wavenumber qu=1.5 for 
T/n=0 (dashed), 0.1 (full line), 2.01 (dash-dotted), and 40.0 (long-dashed). 
Horizontal lines indicate function range, Eq. (l20|) . for this wavenumber. 



this frequency gap has the following meaning. Long-wavelength excitations 
of the trapped ideal gas do not constitute propagating density fluctuations 
(phonons). Instead, a long-wavelength perturbation of the harmonically 
trapped gas will induce a collective oscillation with frequency Q of the total 
gas mass. [12] The appearance of an excitation with frequency Ug could be 
thought of as a hybridization of thermal free-particle excitations and the 
collective oscillation with frequency Q associated with the trap. With in- 
creasing parameter qu, however, this mode will disappear merging into the 
continuous spectrum. 



which coincides with the well-known relaxation spectrum of density fluctu- 
ations of a uniform ideal gas for S{q) = l.[T3] As qu is increased above one, 
the characteristic frequency will change its physical significance, because a 
rapidly growing number of additional resonances at ^^=0, ±2Q, ±3Q, . . . will 
contribute to the second-frequency moment (see Eq. fl28|) below). Eventually, 
for wavenumbers qu ^ 1, the characteristic frequency ujg — > qv^ will deter- 
mine the thermal width of the broad Gaussian spectral line centered at zero 
frequency. Note that the normalized dynamic structure factor has become 




(23) 
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Figure 3: Normalized dynamic structure factors at T=5To. For r/Q=0.1, at 
qu=1.5 (full line) from Eqs.([5]),([23) and qu=0 (dash-dotted) from Eq.([21]). 
Enlarged by factor 8, for r/Q=0, and at qu=1.5: (5-resonances (vertical 
grid lines) and weights (dashed stairway) from Eqs. (!28|) . fl22l) . and large-(gM) 
asymptote (long-dashed) from Eq. (l23l) . 



independent of the trap frequency in the free-particle limit described by 
Eq.m- 

Both S{q,uj) and (j)"iq,uj), Eq.([22]), may be computed from Eqs.(l5|),(l71).(fT8l) 
for arbitrary g > and F > 0. An example is given in Fig. [3]showing 0"(g, u) 
(full line) at wavenumber qu=1.5 of a trapped low-friction gas (F=0.1 Q). It 
corresponds to the intermediate scattering function which is displayed in 
Fig. [2] by a full line, too. Besides the friction-broadened oscillation peak at 
u!=Q, the spectrum clearly exhibits the above mentioned additional peaks 
appearing at uj=0 and 2Q. These additional peaks are absent in the cor- 
responding spectrum taken at qu=0, which is also displayed in Fig. [3l In 
the long-wavelength limit, the normalized dynamic structure factor may be 
calculated analytically, since 

]im 6"{q,u;) ^ 0"(^) := - / dt e'*- 0(t) = . ^ , , (24) 

where 0"(a;) denotes the relaxation spectrum of the oscillator elongation, 
which reduces to a pair of 5-functions for vanishing friction constant: 

\iin(j)"{0,uj) = ^[6{u-Q)+6{uj + Q))] . (25) 
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w/Q 



0.5 1 1.5 2 2.5 3 3.5 



Figure 4: Normalized dynamic structure factors for r=0.1f2, q=1.5uQ , 
and T/To=20, 7 and 5 (dash-dotted, dashed, and full line, resp.) from 
Eqs.(l5]),([22D, and for T/To=5 (long-dashed) from Eq.(l23|). 



For sufficiently small friction, ^"(O,^^) describes the collective density excita- 
tion at frequency Uq=o=Q exhibited in Fig. [3] (dash-dotted curve). It should 
be noted, however, that its detection in a scattering experiment will be dif- 
ficult due to the small overall scattering intensity S{q) ~ {quY associated 
with small values of the parameter qu. On the other hand, one has a large 
overall scattering intensity of S{q) ~ 0.9 for qu=1.5, which corresponds to 
the full curve in Fig. [3] and to all curves in Fig. HI 

At non-zero q and sufficiently small F, one may speculate the spectrum 
(j)"{q,uj) of the form Eq. flMl) but with — > fig and F —>■ Tg, where Qq ~ Uq 
as g ^ 0. 

Obviously, the damping mechanism introduced above cannot create positive 
dispersion {Qq > Q). Weak frictional and stochastic forces merely achieve 
broadening of the resonances at ±fi without shifting them. This broadening 
will only in the over-damped regime lead to an apparent negative dispersion 
associated with the collapse of peaks at u!=0 (see Eq. (Bl|) ). 

Figure m displays variations of the normalized S{q,uj) induced by changing 
the temperature. If uq=\J k-olQ/ (mfi^) serves as unit of length, where Tq 
denotes a reference temperature, all T-dependence of the above formulas will 
be contained in qu=quQjT /Tq. Choosing a fixed wavenumber q=l.buQ^ and 
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the BEC critical temperature of M=2 x 10^ ^^Rb atoms in a trap with 0=27r x 
120 Hz as a reference, i.e. Tq=Tc ~ 147 nK, the full curve in Fig. H] will 
correspond to T=5 Tc ~ 735 nK, a temperature well above the condensation 
point, where the linear extension of the cloud is L 28 fim. 

We note that, for a gas at an ultra-low temperature well above Tc, the 
difference between density-fluctuation spectra of a trapped and a free ideal 
gas is not only qualitative but quantitative as the comparison of full and 
long-dashed curves in Fig. IHmay demonstrate. This makes us believe that 
the predicted spectra could be observed in experiments. 

Although we chose a small damping constant r=0.1f2 in producing FigsjS] 
andlH the friction may be unrealistically large, if the results are interpreted 
as describing an ultra-cold gas cloud. Instead of producing similar plots 
from Eqs.flSl). fl2^ for much smaller friction constant, which would introduce 
an increasing amount of difficulties in handling the numerical integration, we 
will present the analytical solution of the frictionless case in the following 
section. 

5 Undamped motion and DWF 

For a purely harmonic trap, one has (j){t)=cos{Qt) according to Eq. flA.4p . 
Using this in Eq. ffTOl) . one obtains the intermediate scattering function 



where the generating function of modified Bessel functions of the first kind 
[m Chap. 9.6.34 ] was applied in the last line. Noting /_„(z)=/„(2;), the 
dynamic structure factor resulting from Eq. fl271) reads explicitly: 



S{q, u) = Zo{q) 6{uj) + ^ [6{u - nQ) + 6{u + nQ)] . (28) 




(26) 



(27) 



oo 



n=l 



with spectral weights 



Zn{q) : 




Z-n{q) 



(29) 
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0.25r 




Figure 5: Spectral weights of S{q,u!) for n=0 (full), 1 (dashed), 2 (dash- 
dotted), 3 (long-dashed), 4 (long-dash-dotted). 



corresponding to excitation peaks situated at uj=nfl (n=0, ±1, ±2, . . .), which 
sum up to E^=_oo Zn{q) = S{q). 

Weights are displayed in Fig. [5] as function of the dimensionless parameter 
qu. The weights of the elastic and the first inelastic peak at uj= ± Q have 
pronounced maxima, while all other Zn{q) also show a maximum of smaller 
magnitude around qu=2, decrease to zero slowly for g ^ oo, and vanish at 
q=0. For \n\ > 1, they decrease with increasing n. Further, it can be seen 
that Zi clearly dominates the inelastic-scattering weights Zn in the range of 
small wavenumbers, 1 > qu > 0. In the long- wavelength limit, finally, the 
collection of infinitely many resonances will collapse into a pair at uj= ± fl, 
in accordance with Eg. (1251) . That is easily verified by applying the limiting 
values 

Znjq) ^ 1 
9-0 S{q) 2 

in Eg. (1251) from which one finds 



lim^^^ = i5H,i, (30) 



5(g, u) — > gV/2 [5{iu - fi) + 6{iu + n)] (31) 

for qu <^ 1. Spectral weights nZn{q)/ S{q) corresponding to (f)"{q,Lj) have 
been indicated for qu=1.5 in Fig. [3] in terms of a 'stairway' spectral function 
together with the continuous spectrum Eg. (l23|) . Plotting a seguence of such 



13 



f (q) 




Figure 6: Debye-Waller factor of Af=2. x 10^ atoms in harmonic trap for 
r > (full line) and r=0 (dashed). DWF of single trapped atom {J\f=l) for 
r > (long-dashed) and r=0 (dash-dotted). 



stairways for increasing values of the parameter qu will provide a good vi- 
sualization of the spectrum metamorphosis from a sum of only a few visible 
stairs for qu < 1 to the continuous broad spectrum of the uniform ideal gas 
for qu ^ 1. 

It is to be noted that for the purely harmonic trap the dynamical structure 
factor S{q, u) defined by the temporal Fourier transform of the correlation 
function of density fluctuations, Eq.dS]), will contain an elastic scattering con- 
tribution due to F{q, oo) ^ 0, which is absent for F > 0. This contribution, 
together with the a-priori elastic contribution 6{uj) \ (iVq) P/A/" to the coher- 
ent double-differential scattering cross section, when used in Eq.([6]), will lead 
to DWF of the trapped ideal gas 

f(q) - ^-^^"^[^o(^^^^)+^-l] (32) 



whereas for F > the bracket in the numerator of Eq. (1321) is simply replaced 
with A/". Contrary to a homogeneous gas, the DWF of the trapped atomic 
cloud does not vanish and depends on A/". The DWF is displayed for a 
trapped cloud of A/'=20, 000 atoms in Figj6l For comparison, also the DWF 
of a single trapped atom (A/'=l) is given, which coincides with the Lamb- 
MoBbauer factor (LMF) of the trapped ideal gas. 
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6 Concluding remarks 



In this paper we have calculated the density-fluctuation spectrum of a clas- 
sical ideal gas in a harmonic potential at temperatures well above the BEC 
temperature, but in a micro Kelvin range which is still of interest. In this 
regime, HQ is much less than k-sT for trapped atomic gases. It is found that 
the difference between the density fluctuation spectra of the trapped and 
the free gas is large enough to be observed experimentally. The calculation 
predicts the existence of density excitations of the trapped mass which in- 
creases with increasing wave number. It is shown that Debye-Waller factor 
for trapped thermal clouds of J\f atoms are finite contrary to a homogeneous 
gas, opening the possibility of observing recoilless emission of ultra-low fre- 
quency radiation from the trapped gas. It is hoped that this work will stim- 
ulate interest in performing inelastic scattering experiments. 
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Appendix 

A Time evolution of r(^) and p(t) 

In order to evaluate Eq. (ll7p . we need to know the time evolution of particle 
position r(t)=R(p, f, t) and momentum p(t)=P(p, f , t). In the present case 
of an isotropic oscillator potential, the vector components of r{t) as well as 
of p{t) are not coupled. Therefore we outline the method of evaluation of 
the x-components x{t) and Px(t), only. The same method apphes to the y- 
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and 2;-components. 

In the presence of a driving force J-'x{t), the elongation x{t) will obey the 
damped driven oscillator-equations of motion, 

^xit) = -p,it) (A.l) 
dr m 

= -mn^ x{t) - mT ^x(t) + J^.,{t) , (A.2) 
with the unique solution 

x{t) = (Pit -to)x + ixit - to) Px + f'dt' ixit - t') J^x{t') (A.3) 

corresponding to initial conditions x{tQ)=x and Px(to)=px. Here (j){t) and 
X(t), defined by {t > 0) 

0(t) := ^ ^ , X(t):=: , (A.4) 

Zi - Z2 im Zi - Z2 



zi,2 = -r/2±^(r/2)2-fi2^ (A.5) 

denote the (normalized) relaxation function and the linear response function 
of the oscillator elongation, respectively. Equations similar to Eq. (lA.3p also 
apply for the y- and 2;-components. 

Indeed, (j){t) represents the normalized Kubo relaxation function of elonga- 
tion. This can be seen as follows. Choosing to=0 in Eq. flA.3p . multiplying 
with x/{kBT) on both sides and subsequently taking a thermal average will 
— on the left-hand side — result in the Kubo relaxation function of elon- 
gation, $(t)=(x(t)x) /{k-QT). On the right-hand side, only the first term 
(pit) (x^) /(/cbT) = 0(t)$(O) will survive, where $(0)=l/(mfi2). 

Both functions, 0()f:) and x(t), are related by Kubo's identity, 

which can be verified explicitly for the results given in Eq.dAH). Also, the 
initial values of the relaxation function, 

0(0) = 1 , 0(0) = , (A.7) 
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can be read from Eqs. (lA.4p . (lA.6p . while one has 



oo 







oo 







(r>o) 



(A.8) 



for finite friction. The magnitude of the normahzed relaxation function is 
restricted to the interval [0, 1] or, equivalently, 



- 1 < < 1 



(A.9) 



for any time t. 

In addition to the restoring force and the friction force, it has been assumed 
in Eq.( ]A.2l) that there acts a time-dependent random force J-'x{t) on the oscil- 
lator. Thermal averaging is then understood to include additional averaging 
over stochastic variables as indicated explicitly in Eq. (fT7I) . For evaluating 
the stochastic average, we assumed that the driving force J-'x{t) represents 
Gaussian white noise, which after straightforward calculation leads to the 
result Eq.([T5|). 



B Non— diagonal correlation functions 

Assuming undamped oscillations 

- / N X ^ sin(i7t) „ 

f (t) = cos(fit) f + ]^ ' p 



one finds from Eqs.( fT6] - [T7ll 

F„,(t)=e 



mil 



which by Fourier transformation leads to 
F(r,r',t) : 



I fr-r'cos(m)]M / ^2 



(m| sin(l]t)|v/27r) (u^/2^T 



exp 



u\/2n 



(B.l) 



(B.2) 



(B.3) 



It is noted that Fqq/(t)=Fq/q(t) and limg^o -^qq'(^)=0- W^e further observe 
/dV F(r, r', t) = /dV F(r, r', t)=0. 



17 



Taking the limit t ^ in Eq. (lB.3p one finds 



F(r,r',0) 



o(r — r ) 



(B.4) 



which is the inverse Fourier transform of Eq. (fT2!) or of the it=0-value of 
Eq.dRl. 
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